Suzaku Observation of Group of Galaxies NGC 507: 
Temperature and Metal Distributions 
in the Intra-cluster Medium 
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Abstract 

Temperature and abundance distributions of the intra-cluster medium (ICM) in the NGC 507 group 
of galaxies were studied with Suzaku. Observed concentric annular spectra were well-represented by a two 
temperature model for ICM, and we found steeper abundance gradients for Mg, Si, S, and Fe compared 
with O in the central region. Abundance ratios of ev-elements to iron were found to be similar to those in 
other groups and poor clusters. We calculated metal mass-to-light ratios for Fe, O and Mg (IMLR, OMLR, 
MMLR) for NGC 507, and values for different systems were compared. Hotter and richer systems tend to 
show higher values of IMLR, OMLR, and MMLR. OMLR and MMLR were measured to an outer region 
for the first time with Suzaku, while IMLR was consistent with that with ASCA. We also looked into 
2-dimensional map of the hardness ratio, but found no significant deviation from the circular symmetry. 

Key words: galaxies: clusters: individual (NGC 507), galaxies: intergalactic medium, galaxies: 
abundances 
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1. Introduction 

Recent X-ray observations have shown metal abundance 
profiles in the ICM based on the spatially resolved spec- 
tra. Groups of galaxies carry important information about 
the entire history of cosmic chemical evolution, since the 
groups contain 50-70% of galaxies in the universe (e.g. 
Mulchaey et al. 2006) and they serve as building blocks 
of clusters of galaxies. In order to know how the ICM 
has been enriched, we need to measure the amount and 
distribution of metals in the ICM. Because Si and Fe are 
both synthesized in type la and type II supernovae (SNe 
la and II), we need to know O and Mg abundances, which 
are synthesized practically in SNe II, in resolving the past 
metal enrichment process in ICM by supernovae in early- 
type galaxies (Arnaud et al. 1992; Renzini et al. 1993). 
Comparison of the ICM properties between groups and 
clusters of galaxies will also enable us to understand the 
bottom-up process of the hierarchical structure formation 
in the universe. 

ASCA measured the distribution of the heavy elements, 
such as Si and Fe, in the ICM (Fukazawa et al. 1998; 
Fukazawa et al. 2000; Finoguenov et al. 2000; Finoguenov 
et al. 2001). Renzini (1997) and Makishima et al. (2001) 
summarized iron-mass-to-light ratios (IMLR) for various 
objects, as a function of their plasma temperature serv- 
ing as a measure of the system richness with ASCA, and 



they showed that the early-type galaxies released large 
amount of metals formed through past supernovae explo- 
sions as shown earlier by Arnaud et al. (1992). Recent 
XMM-Newton and Chandra observations have enabled us 
to study properties of the heavy elements in the ICM in 
detail. These observations showed not only Si and Fe, 
but also O and Mg abundance profiles, however the O 
and Mg abundance measurements have been limited only 
for the central regions of very bright clusters or groups 
of galaxies dominated by cD galaxies in a reliable man- 
ner (Finoguenov et al. 2002; Xu et al. 2002; Matsushita 
et al. 2003; Tamura et al. 2003; Buote et al. 2003a; Buote 
et al. 2003b; Humphrey & Buote 2006). Tamura et al. 
(2004) derived IMLR for five clusters, and the oxygen 
mass for several clusters with XMM-Newton. However, 
oxygen-mass-to-light ratios (OMLR) for rich clusters are 
not reliable due to the lower emissivity of O VII and O VIII 
lines in high temperatures. Matsushita et al. (2007b) 
showed different distribution profiles among the metals, 
which suggest difference in the process of metal enrich- 
ment. The abundance measurements of O and Mg with 
XMM-Newton, particularly for the outer regions of groups 
and clusters, are quite difficult due to the relatively high 
intrinsic background. Suzaku XIS can measure all the 
main elements from O to Fe, because it realizes lower back- 
ground and higher spectral sensitivity, especially below 1 
keV (Koyama et al. 2007). Matsushita et al. (2007a), Sato 
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ct al. (2007a), Sato et al. (2008), and Tokoi et al. (2008) 
have shown the abundance profiles of O, Mg, Si, S, and 
Fe with Suzaku to the outer regions with good accuracy. 

NGC 507 is a nearby group of galaxies (z = 0.01646) 
characterized by a smooth distribution of ICM. The 
ICM properties have been studied with ROSAT (Kim & 
Fabbiano 1995; Paolillo et al. 2003), ASCA (Matsumoto et 
al. 1997), Chandra (Kraft et al. 2004; Humphrey & Buote 
2006; Rasmussen & Ponman 2007), and XMM-Newton 
(Kim & Fabbiano 2004). Kim & Fabbiano (1995) revealed 
a cooler central region with ROSAT PSPC, and Kim & 
Fabbiano (2004) showed the supersolar metal abundances 
within the D 2 5 ellipse of NGC 507 with XMM-Newton. 
On the other hand, Humphrey & Buote (2006) showed a 
near solar metal abundances, the same as derived by Kraft 
et al. (2004) with Chandra. Kim & Fabbiano (2004) and 
Humphrey & Buote (2006) showed the a-elements to iron 
ratios, and both the abundance ratios were almost consis- 
tent to be ~ 1 solar. Kraft et al. (2004) also reported a 
sharp edge or a discontinuity in the radial surface bright- 
ness profile 55 kpc east and southeast of NGC 507 covering 
an ~ 125° arc. NGC 507 is also a known radio source and 
has been classified as an FR I radio galaxy (Fanti et al. 
1986), and the edge of the radio lobe corresponds to the 
discontinuity of the X-ray surface brightness as shown in 
Kraft et al. (2004). 

This paper reports on results from Suzaku observa- 
tions of NGC 507 out to 13' ~ 260 kpc, correspond- 
ing to ~ 0.24 riso . We use H = 70 km s _1 Mpc~\ 
= 1 — Qjvf = 0.73 in this paper. At a redshift of 
z = 0.01646, 1' correspo nds to 20.1 k pc, and the virial 
radius, r 180 = 1.95 h^ ^k{T)/ 10 keV Mpc (Markevitch 
et al. 1998), is 1.08 Mpc (54') for an average temper- 
ature of k(T) = 1.5 keV. Throughout this paper we 
adopt the Galactic hydrogen column density of Ah = 
5.24 x 10 20 cm~ 2 (Dickey & Lockman 1990) in the di- 
rection of NGC 507. Unless noted otherwise, the solar 
abundance table is given by Anders & Grevesse (1989), 
and the errors are in the 90% confidence region for a sin- 
gle interesting parameter. 

2. Observations and Data Reduction 

2.1. Observation 

Suzaku observed the central region of NGC 507 in July 
2006 (PI: K. Sato) with an exposure time of 79.6 ks. The 
observation log is given in table 1, and the XIS image in 

0. 5-4 keV is shown in figure 1. We analyze only the XIS 
data in this paper, although Suzaku observed the object 
with both XIS and HXD, which acquired the data. The 
XIS instrument consists of four sets of X-ray CCDs (XIS 0, 

1, 2, and 3). XIS 1 is a back-illuminated (BI) sensor, while 
XIS 0,2, and 3 are front-illuminated (FI). The instrument 
was operated in the Normal clocking mode (8 s exposure 
per frame) , with the standard 5 x 5 or 3 x 3 editing mode. 

It is known that the optical blocking filters (OBF) of the 
XIS have gradually been contaminated by outgassing from 
the satellite. The thickness of the contaminant is different 
among the sensors, and is also dependent on the location 
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Fig. 1. Combined XIS image of central observation in 
the 0.5—4.0 keV energy range. The observed XISO-3 im- 
ages were added on the sky coordinate after removing 
each calibration source region, and smoothed with <r = 16 
pixel ~ 17" Gaussian. Estimated components of ex- 
tragalactic X-ray background (CXB) and instrumental 
background (NXB) were subtracted, and the exposure 
was corrected, though vignetting was not corrected. 



Table 2. Estimated column density of the contaminant for 
each sensor at the center of CCD in units of 10 18 cm -2 . 





XIS0 


XIS1 


XIS2 


XIS3 


Carbon 


2.52 


3.94 


3.83 


5.78 


Oxygen 


. ... 0.419 


0.657 


0.638 


0.963 



on the CCD chips in the way that it is thickest in the 
center. Estimated column densities (C/0=6 in number 
ratio is assumed) during the observation at the center of 
the CCD are given in tabic 2. 1 We included these effects in 
the calculation of the the Ancillary Response File (ARF) 
by the "xissimarfgen" ftools task of 2006-10-26 version 
(Ishisaki et al. 2007). Since the energy resolution also 
slowly degraded after the launch, due to radiation damage, 
this effect was included in the Redistribution Matrix File 
(RMF) by the "xisrmfgen" Ftools task of the 2006-10-26 
version. 

2.2. Data Reduction 

We used version 1.2 processing data (Mitsuda et al. 
2007), and the analysis was performed with HEAsoft ver- 
sion 6.1.1 and XSPEC 11.3.2t. The analysis method was 
almost the same as those in Sato et al. (2007a) and Sato 
et al. (2008). However, because the observations were 
not supported by the Good-Time Intervals (GTI) defined 
to exclude the telemetry saturation by the XIS team, we 

Calibration database file of ae_xiALcontami_20061024 . f its was 

used for the estimation of the XIS contamination (N = 0, 1,2,3 
denotes the XIS sensor). 



Suzaku observation of NGC 507 

Table 1. Suzaku Observation logs for NGC 507. 
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Object 


Seq. No. 


Obs. date 


(RA, Dec)* 


Exp. 


After screening 








J2000 


ksec 


(BI/FI) ksec 


NGC 507 


801017010 


2006-07-28T14:51:16 


(01 h 23 m 40?0, +33°15'21") 


79.6 


79.2/79.2 



* Average pointing direction of the XIS, written in the RA_NOM and DECJMOM keywords of the event FITS files. 



could not execute the GTI correction. The light curve 
of each sensor in the 0.3-10 keV range with a 16 s time 
bin was also examined in order to exclude periods with 
anomalous event rates which were greater or less than 
±3ct around the mean. After the above screenings, the 
remaining exposure time of the observation stayed almost 
unchanged as shown in table 1. The exposure after the 
screening was essentially the same as that before screening 
in table 1, which indicates that the non X-ray background 
(NXB) was almost stable during the observation. Event 
screening with cut-off rigidity (COR) was not performed 
in our data. 

In order to subtract the NXB and the extra-galactic 
cosmic X-ray background (CXB), we employed the dark 
Earth database of 770 ks exposure, provided by the XIS 
team for the NXB, and employed the CXB spectrum 
given by Kushino et al. (2002). These analysis meth- 
ods were also the same as in Sato et al. (2007a) and 
Sato et al. (2008). In order to remove influence from 
bright sources, we eliminated the regions within 80 arcsec 
from (01 h 23 n T4?0, +33°10'57") and within 45 arcsec from 
(01 h 23 m 39fl, +33 o 09'04") in our field of view. However, 
because the spatial resolution of Suzaku is not enough to 
resolve all point sources especially at the central region, 
effect from the remaining weaker point sources was taken 
into account in the spectral fitting function in the next 
sections. 

2.3. Generation of ARFs 

A precise surface brightness profile of NGC 507 was 
needed to generate the Suzaku ARF, and we used the 
XMM-Newton image which had much better spatial reso- 
lution than Suzaku. For this XMM-Newton data, we used 
MOS1+2 data (34.7 ks), and followed the data reduction 
by Sato et al. (2005). We subtracted blank-sky data as 
the background (Read & Ponman 2003), eliminated point 
sources in the ICM, and corrected the vignetting. Figure 2 
shows a radial profile of NGC 507 with XMM-Newton in 
the 0.8-3 keV range, fitted with a j3 model. The origin of 
the profile is set at (RA, Dec) = (l h 23 m 39?7, +33°15'22") 
in J2000.0. The best-fit parameters are j3 = 0.48 and r c = 
l'.Q. Though the fit is not acceptable, the ARF showed 
little influence on the temperature and abundance in the 
spectral fits. We then generated two ARFs for the spec- 
trum of each annular sky region, A and A B , which respec- 
tively assumed uniform sky emission and ~ 0.5° x 0.5° size 
of the /3-model surface brightness profile obtained with 
the XMM-Newton data. We did not use the raw XMM- 
Newton image, but the smoothed image derived from the 
parameters of the /3-model fit to generate the ARFs, be- 
cause the raw image had a gap between the CCD chips. 
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Fig. 2. In the upper panel, a radial profile of the surface 
brightness of NGC 507 in the 0.8-3 keV band is plotted for 
XMM-Newton MOS1+2 (r < 10'). The best-fit model is 
shown by the solid gray line. In the bottom panel, the fit 
residuals are shown in units of cr, which correspond to the 
data minus the folded model, divided by la error of each 
data point. 

In this way, NGC 507 was characterized by a smooth and 
symmetric ICM distribution. 

3. Temperature and Abundance Profiles 

3.1. Spectral Fit 

We extracted spectra from five annular regions of 0'- 
2', 2'-4', 4'-6', 6'-9', r > 9', centered on (RA, Dec) = 
(l h 23 m 40?0, +33°15'21"). Table 3 lists the areas of the 
extracted regions (arcmin 2 ), fractional coverage of the an- 
nulus (%), the SOURCE_ratio_REG values (%; see caption 
for its definition) and the BI and FI counts for the ob- 
served spectra and the estimated NXB and CXB spec- 
tra. The fraction of the background, /bgd = (NXB + 
CXB)/OBS, was less than ^30% even at the outermost 
annulus, although the Galactic component is not consid- 
ered here. Each annular spectrum is shown in figure 3. 
The ionized Mg, Si, S, Fe lines are clearly seen in each 
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Table 3. Area, coverage of whole annulus, SOURCE_RATIO_REG, and observed/estimated counts for each annular region. 



Region 


* Area T 


Coverage 


T source.* Energy 




BI counts § 






FI counts & 






(arcmin 2 




RATIO.REG (keV) 


OBS 


NXB 


CXB 


/bgd 


OBS 


NXB 


CXB 


/bgd 


0'-2' 


12.6 


100.0% 


29.2% 0.4-7.1 


23,452 


335 


338 


2.8% 


49,445 


509 


794 


2.6% 


2'-4' 


37.7 


100.0% 


21.9% 0.4-7.1 


26,997 


1,023 


957 


7.3% 


53,622 


1,501 


2,284 


7.1% 


4'-6' 


62.8 


100.0% 


14.4% 0.4-7.1 


18,274 


1,708 


1,405 


17.0% 


34,960 


2,488 


3,356 


16.7% 


6'-9' 


133.3 


94.3% 


13.4% 0.4-7.1 


21,286 


3,490 


2,419 


27.8% 


38,891 


4,961 


5,520 


26.9% 


r >9' 


64.0 


23.1% 


4.0% 0.4-3.0 


5,481 


927 


737 


30.4% 


9,420 


1,218 


1,497 


28.8% 



SOURCE_RATIO_REG represents the flux ratio in the assumed spatial distribution on the sky (/3-modcl) inside the accumulation region to 
the entire model, and written in the header keyword of the calculated ARF response by "xissimarfgen" . 

* The outermost region (r > 9') include the calibration source region, 
t The largest values among four sensors are presented. 

* SOURCE.RATIO.REG = COVERAGE X J^° ut S(r) rdr/ S(r) rdr, where S(r) represents the assumed radial profile of NGC 507, and 
we defined S(r) in 30' X 30' region on the sky. 

§ OBS denotes the observed counts including NXB and CXB in 0.4-7.1 keV or 0.4-3 keV. NXB and CXB are the estimated counts. 




0.1 I . 1 , 1 , 1 1 

0.2 0.4 0.6 

Oxygen abundance of ICM (solar) 

Fig. 4. A plot of confidence contour between kT of apec 
component (i.e. Galactic component) and the O abun- 
dance of vapec model (i.e. ICM) for the r > 9' annu- 
lus, in the simultaneous fitting of all annuli with the 
apec + phabs X (vapec 1 + vapec 2 + zbremss) model. The 
cross denotes the best-fit location, and the two contours 
represent 1<t and 90% confidence ranges, from inner to 
outer, respectively. 

region. The OVII and OVIII lines are prominent in the 
outer rings, however, most of the O VII line is considered 
to come from the local Galactic emission, and we dealt 
with those in the same way as Sato et al. (2007a) and 
Sato et al. (2008). 

The spectra from BI and FI for all regions were fitted si- 
multaneously. For the four inner regions the fitted energy 
range was 0.4-7.1 keV, and it was 0.4-3.0 keV for the 
outermost region. In the simultaneous fit, the Galactic 
emission component as a background was common for all 
regions, while the NGC 507 emission component was un- 
linked between each region. We excluded the narrow en- 
ergy band around the Si K-edge (1.825-1.840 keV) be- 
cause of incomplete response. The energy range below 
0.4 keV was also excluded because the C edge (0.284 keV) 
seen in the BI spectra could not be reproduced well in our 
data. The range above 7.1 keV was also ignored because 
Ni line (~ 7.5 keV) in the background left an artificial 
structure after the NXB subtraction at large radii. In the 



simultaneous fit of BI and FI data, only normalization 
parameters were allowed to be different between them, 
although we found that the derived normalizations were 
quite consistent to be the same. For the outermost region 
(r > 9') we fitted the spectrum in 0.4-3.0 keV to exclude 
contribution from the Fe-55 calibration source. 

It is important to estimate the Galactic component pre- 
cisely. In order to determine surface brightness and spec- 
tral shape of the Galactic component, we carried out si- 
multaneous fit for all annuli. The Galactic component 
gives significant contribution in these annuli as shown in 
figure 3, however the ICM component is still dominant in 
almost all the energy range except for the OVII line. 

We assumed either one or two temperature apec model 
for the Galactic component, and fitted the data with 
the following model formula: apec 1T or 2T + phabs x 
(vapec 1T or 2T 4- zbremss). The resultant normalizations 
of the apec models in table 4 are scaled so that they give 
the surface brightness in the unit solid angle of arcmin 2 , 
and the normalizations of the apec components were con- 
strained to give the same surface brightness and the same 
temperature for all annuli. When we fitted the data with 
two apec models whose temperatures were free parame- 
ters, we could not obtain reasonable fit. Thus, we fixed the 
two temperatures to be 0.1 and 0.3 keV, which were con- 
sistent with the previous studies (e.g. Lumb et al. 2002). 
When we assumed two apec model, whose temperatures 
were fixed to be 0.1 and 0.3 keV, the resultant parameters 
were almost the same as those in table 5. As a result, 
we concluded that only one apec model was enough to fit 
the NGC 507 data since the fit improvement showed low 
significance with the two apec model. Therefore, we per- 
formed simultaneous fit in the 0.4-7.1 keV range (exclud- 
ing 1.825-1.840 keV) for the four inner regions and in the 
0.4-3.0 keV range for the r > 9' region, assuming an apec 
model for the Galactic component, by following model, 
apec + phabs x (vapec 1 + vapec 2 + zbremss), as shown in 
table 4. 

To demonstrate how sensitive the O abundance of the 
ICM is to the assumed Galactic component model, we 
show in figure 4 a confidence contour between kT (keV) 
of the apec component and the O abundance (solar) of 
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Fig. 3. The panels show the observed spectra at the annular regions of NGC 507 which are denoted in the panels, and they are 
plotted by red and black crosses for BI and FI, respectively. The estimated CXB and NXB components are subtracted, and they are 
fitted with the apec + phabs X (vapec 1 + vapec 2 + zbremss) model shown by green and yellow lines for the BI and FI spectra. The 
vapec 1 (hot) and vapec 2 (cool) components of ICM correspond to blue and magenta lines, respectively. The apec 
component for the BI spectra are indicated by cyan line. The zbremss component for the BI spectra within r < 4' 
region are indicated by orange line. The energy range around the Si K-edge (1.825—1.840 keV) is ignored in the 
spectral fit. The lower panels show the fit residuals in units of a. 



Table 4. The best-fit parameters of the apec component for the simultaneous fit of all spectra of NGC 507 with one or two 
temperature models (apec) for Galactic emissions and phabs X (vapec 1T or 2T + zbremss r<i i) model for ICM. 



Fit model Norrm * Wi Norm 2 * kf 2 x 2 /dof 
(keV) 



apec 1 - 


- phabs 


x (vapec 1 


+ zbremss) 


0.47 


0.110 






3173/1520 


apec 1 - 


- phabs 


x (vapec 1 


+ vapec 2 + zbremss) 


0.42 ±0.14 


n 1 KQ+0.034 






2067/1510 


apec 1 - 


- apec 2 


+ phabs x 


(vapec 1 + zbremss) 


0.26 


0.1 (fix) 


0.33 


0.3 (fix) 


3073/1520 


apec 1 - 


- apec 2 


+ phabs x 


(vapec 1 + vapec 2 + zbremss) 


1.02 ±0.42 


0.1 (fix) 


0.13 ±0.13 


0.3 (fix) 


2075/1510 



* Normalization of the apec component divided by the solid angle, fi u , assumed in the uniform-sky ARF calculation (20' radius), 
Norm = / n a n n dV I (4tt (1 + z) 2 D 2 ) / Q' [ xlO 20 cm 5 arcmin 2 , where Da is the angular distance to the source. 
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Fig. 5. (a): Radial temperature profiles derived from the spectral fit for each annulus against the projected radius. Black dashed 
lines show systematic change of the best-fit values by varying the thickness of the OBF contaminant by ±10%. Light-gray dashed 
lines denote those when the estimated CXB and NXB levels are varied by ±10%. (b)-(f): Radial abundance profiles derived and 
plotted in the same way as in (a), (f): XMM results of Fit 13 in Kim & Fabbiano (2004) correspond to gray crosses. 



vapec for the outermost annulus (r > 9'). There seems a 
negative correlation between the two parameters, because 
higher temperature of the Galactic component produces 
more intense O VIII line relative to O VII one, which pushes 
to reduce the O VIII line from the ICM (vapec component). 
Influence on the derived temperature and abundance by 
the modeling of the Galactic component will be examined 
in subsection 3.2, too. 

The ICM spectra for all regions were clearly better rep- 
resented by two vapec models than one vapec model as 
shown in table 5. The abundances were linked in the 
following way, Mg=Al, S=Ar=Ca, Fc=Ni, and were also 
linked between the two vapec components for each region. 
The fit results are shown in table 5. Kim & Fabbiano 
(2004) showed that the model for the central region needed 
to include the low-mass X-ray binary (LMXB) component 
based on the ROSAT and XMM-Newton observations, we 
included a kT — 7 keV zbremss model for this component 
in the fit for 0-2' and 2-4' annuli. The resultant flux in 
0.3-8 keV was - 1 x 10~ 13 and - 4 x 10~ 14 erg s _1 cur 2 
for 0-2' and 2-4' region, respectively, and the flux due to 
the LMXB component with Suzaku was consistent with 
the Kim & Fabbiano (2004) result. 

Although the fit was not statistically acceptable due 
mainly to the very high photon statistics compared with 
the systematic errors in the instrumental response, the 
results were useful to assess whether each element abun- 
dance was reasonably determined or not. Abundance of 
Ne was not reliably determined due to the strong and com- 
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Fig. 6. The ratio of emission measure of the hot ICM 
component to that of the cool ICM component. 

plex Fe-L line emissions, however we left these abundance 
to vary freely during the spectral fit. 

Results of the spectral fit for each annulus arc summa- 
rized in table 6 and figure 4, in which systematic error 
due to the OBF contamination and background (CXB + 
NXB) estimation are included. We examined the results 
by changing the background normalization by ±10%, and 
the error range is plotted with light-gray dashed lines in 
figure 4. The systematic error due to the background es- 
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Table 5. Summary of the parameters from one or two vapec fit (apec + phabs X (vapec 1T or 2 t + zbremss r< 4i )) to each annular 
spectrum of NGC 507. All annuli were simultaneously fitted. Errors are 90% confidence range of statistical errors, and do not 
include systematic errors. The solar abundance ratio of angr was assumed. These results are plotted in figure 5. 



ICM1T 


Normi * 


(keV) 


Norm 2 * 


kT 2 
(keV) 


Normi 1 Norm 2 


x7dof 


0-2' 


527.0 


1.21 


_ 




_ 


1118/339 


2-4' 


99.6 


1.27 


_ 


_ 


_ 


620/339 


4-6' 


29.0 


1.26 


_ 


_ 


_ 


560/340 


6-9' 


9.3 


1.17 


— 


- 


— 


638/340 


r>9' 


3.8 


1.11 


— 


— 


— 


237/162 


total 












3173/1520 


ICM1T 





Ne 


Mg,Al 


Si 


S,Ar,Ca 


Fe,Ni 




(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


0-2' 


0.23 


0.00 


0.25 


0.31 


0.40 


0.22 


2-4' 


0.26 


0.23 


0.33 


0.31 


0.43 


0.26 


4-6' 


0.35 


0.05 


0.56 


0.55 


0.63 


0.42 


6-9' 


0.36 


0.00 


0.59 


0.68 


0.80 


0.51 


r > 9' 


0.21 


0.00 


0.22 


0.19 


0.41 


0.18 


ICM2T 


Normi * 


kTi 
t^Ke v ) 


Norm 2 * 


kT 2 
^Ke v ) 


Normi 1 Norm 2 


X /dof 


0-2' 


115.4 ±7.0 


i - o4 -0.02 


38.0 ±3.2 


o.83±g:8± 


3.04 ±0.32 


576/337 


2-4' 


49.0 ±2.8 


i co+0.04 
i.dO_ 02 


11.3 ± 0.8 


0.85t° ;°i 


4.34 ±0.39 


436/337 


4-6' 


26.1 ±1.7 


1 ci +0.05 
1 - ol -0.04 


5.7±0.6 


0.84l8:g 


4.58 ±0.57 


469/338 


6-9' 


12.6 ±0.9 


i ' DO -0.06 


4.2 ±0.4 


o.siiS:g 


3.00 ±0.36 


399/338 




8 t;xi o 
o . o m l .z 


1 4Q+0-17 


a i -i- n s 

4:. i nz u. o 


n o 9 +0.03 




1 87 /1 fin 

io i j iuu 


total 












2067/1510 


ICM2T 





Ne 


Mg,Al 


Si 


S,Ar,Ca 


Fe,Ni 




(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


0-2' 
2-4' 
4-6' 
6-9' 
r>9' 


o.58±»:i« 

0.50+fe 
36+° 15 
34+ 17 
22+ ' 27 


2 02 +u - 41 

1 34+0-20 

1 oo+°- 28 

±.UU_ 26 
1 1O+0-32 

0.39tg-j§ 


i.osii! 

67+ - 15 
0.58+° ;i 


i oo+O.lS 
1 -°°-Q.13 

OO+0.09 

n 4^+0-os 

u - 4o -0.07 

0-54±g;lg 
0.3418-.S 


1.36+|" 
0.82+°;}° 
^i+ - 10 
^2+ 013 

u - oz -0.13 

^3+°- 31 

U.OO_ 29 


96 +(J U9 
61+ - 05 
40+ - 05 
n 41+0-05 

3^+0- 10 



* Normalization of the vapec component scaled with a factor of S0URCE_RATI0_REG / AREA in table 3, 
Norm= SOURCE.RATIO.REG J ne „ H£i y/[4 7r (i +2 )2 D 2] xl0 -20 cm -S arcm i n -2 ) where Da is thc an g U l ar 

distance to the source. 

t All regions were fitted simultaneously. 



timation is almost negligible. Thc other systematic error 
concerning the uncertainty in the OBF contaminant is 
shown by black dashed lines. A list of x 2 /dof is presented 
in table 6. 

3.2. Temperature Profile 

As mentioned in the previous section, the spectral fits 
needed two temperature model for ICM rather than one 
temperature model. As shown in table 5, the metal abun- 
dances changed dramatically whether we use one or two 
temperature models at the central region. Radial tem- 
perature profile and the ratio of the vapec normalizations 
between the hot and cool ICM components are shown 
in figure 5(a) and figure 6, respectively. Our results for 
the two temperature ICM model are consistent with the 
XMM-Newton results (Kim & Fabbiano 2004) for thc cen- 
tral region within 4'. In the Chandra case, Humphrey & 
Buotc (2006) fitted the spectra with IT model. The pre- 



vious XMM-Newton and Chandra results neglected the 
Galactic component, while our analysis took it into ac- 
count. We can crudely approximate the two temperatures 
as fcTjiot ~ 1-5 keV, and kTc 00 \ ~ 0.8 keV, respectively. 
The cool component is strongest in the innermost region, 
and even though it seems to decline in the outer regions, 
the possible coupling with the Galactic emission makes the 
precise estimation difficult. Thc radius of 10' ~ 200 kpc 
corresponds to <~ 0.19 r i80 , and the temperature decline, 
observed in several other clusters, is not clearly recog- 
nized in this system due partly to the multi-phase nature 
of ICM. The temperatures of cool and hot components 
are consistent with the results of two temperature model 
for NGC 5044 group with XMM-Newton in Buote ct al. 
(2003a). 

Systematic error caused by the background subtrac- 
tion and by the estimation of the XIS filter contamina- 
tion was estimated. We varied the sum of the NXB and 
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Table 6. List of x 2 / dof for each fit of NGC 507. 



Region 


nominal 


contaminaut 


background 






+ 10% -10% 


+10% -10% 


All 


. 2067/1510 


2004/1510 2148/1510 


2104/1510 2094/1510 



CXB flux and the contamination thickness individually 
by ±10%. The results are shown by dashed lines in fig- 
ure 5(a). Though kT coo \ component is not affected by 
these uncertainties, fcTuot shows a significant dependence 
on the background uncertainty most notably in the 6'-9' 
annulus. The tendency looks quite similar to the result 
for HCG 62 with Suzaku (Tokoi et al. 2008). 

3.3. Abundance Profiles 

Metal abundances are determined for the six element 
groups individually as shown in figures 5(b)-(f). O abun- 
dance is strongly affected by the Galactic emission as 
shown figure 4. However, the use of one or two tempera- 
ture model for the Galactic emission did not significantly 
affect the abundance results as mentioned in subsection 
3.1. We also note that Ne abundance has a problem 
in the spectral fit due to the coupling with Fe-L lines. 
Therefore, regarding the spatial structure, we dealt with 
the remaining four elemental groups: Mg, Si, S and Fe. 
The four abundance values and their variation with radius 
look quite similar to each other. The central abundances 
lie between 1.0-1.6 solar, and the abundance decline to 
about 1/5 of the central value is commonly seen in the 
r > 9' annulus. 

Kim & Fabbiano (2004) analyzed XMM-Newton data 
and employed the solar abundance by Grevesse & Sauval 
(1998) giving [Fe/H] = 3.16 x 10~ 5 . This was scaled here 
by a factor of 0.7 to match the Anders & Grevesse (1989) 
value of [Fe/H] = 4.68 x 10~ 5 . As shown in figure 5 (f), 
our result shows lower Fe abundance than the Kim & 
Fabbiano (2004) value within r < 2', while both results 
are almost consistent in 2' < r < 5' region. The spec- 
tral model is the same for both analysis, employing the 
two temperature ICM model with the LMXB component. 
The difference is that our result was derived for the pro- 
jected spectrum while Kim & Fabbiano (2004) used the 
deprojected spectra. The Fe abundance with Suzaku is 
also consistent with the Chandra result by Humphrey 
& Buote (2006). The resultant abundance profiles with 
Suzaku for NGC 507 are also consistent with those with 
the two temperature model for NGC 5044 with XMM- 
Newton in Buote et al. (2003b). 

Again, we looked into the effect of the error by the 
NXB and CXB intensities and the OBF contamination. 
As shown by dashed lines in figures 5(b)-(f), the system- 
atic effect is less than the statistical error for all regions. 

3-4- Direct Comparison of O VII and O VIII Intensities 

We examined the surface brightness of the OVII and 
OVIII emission lines in order to look into the origin of 
the O lines directly from the line intensities. The surface 




Projected radius (') 



Fig. 7. Intensities of O VII and O VIII lines at each an- 
nulus of NGC 507 in units of photons cm~ 2 s _1 sr _1 . 
Numerical values are shown in table 7. 

brightness of the lines was derived by fitting the annular 
spectrum with a power-law + gaussian + gaussian model. 
In this fit, we fixed the Gaussian a to be 0, and allowed 
the energy center of the two Gaussians to vary within 555- 
573 eV or 648-658 eV for O VII or O VIII line, respectively. 
The derived line intensities are summarized in table 7 and 
figure 7. There is a clear excess of the OVIII intensity 
towards the cluster center, while O VII one is consistent to 
be constant. This is a clear evidence that the OVIII line 
is associated with the ICM itself, on the other hand, O VII 
may be due mainly to the Galactic origin. 

4. Discussion 

4--1- Metallicity Distribution in ICM 

The present Suzaku observation of NGC 507 showed 
abundance distribution of O, Mg, Si, S, and Fe out to a 
radius of 13' ~ 260 kpc as shown in figure 5. Ne abundance 
has large ambiguity due to the strong coupling with Fe-L 
lines. Distributions of Mg, Si, S, and Fe are quite similar 
to each other, while O profile in the outer region has a 
large uncertainty. We plotted abundance ratios of O, Mg, 
Si, and S over Fe as a function of the projected radius in 
figure 8. The ratios Mg/Fe, Si/Fc and S/Fe are consistent 
to be a constant value around 1.5-2, while O/Fe ratio for 
the innermost region (r < 2') is significantly lower around 
0.6. In addition, the O/Fe ratio suggests some increase 
with radius. 

Recent Suzaku observations have presented abundance 
profiles in several other systems: an elliptical galaxy 
NGC 720 (Tawara et al. 2008), the Fornax cluster and 
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Table 7. Intensities of O VII and O VIII lines for each annulus of NGC 507 field in units of photons cm -2 s _1 sr — 1 . These intensities 
are derived from the spectral fit with the model power-law + gaussian + gaussian, assuming the uniform-sky ARF response. Results 
by McCammon et al. (2002) measured with a high resolution microcalorimeter array for a large sky area of ~ 1 sr are also presented. 



Region 




Ovu 


vin 


NGC 507 


(0-2') 




26.0±8.4 


NGC 507 


(2-4') 


8.9±4.0 


11.3±4.0 


NGC 507 


(4-6') 


2.5±2.5 


3.1±2.4 


NGC 507 


(6-9') 


4.4±1.7 


3.9±1.6 


NGC 507 


fr>9') 


3.7±2.7 


3.3±2.4 


Galactic average (McCammon et al. 2002) 


4.8±0.8 


1.6±0.4 



0.5 
I 0.2 

£ o.i 

0.05 
1.5 



u 
6! 
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2 
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Fig. 8. Comparison of the abundance ratio for NGC 507 
(red) with results for the HCG 62 group (black; Tokoi et al. 
2008) and the Fornax cluster (gray; Matsushita et al. 2007b). 
Fe abundance and the O/Fe, Mg/Fe, Si/Fe, and S/Fe abun- 
dance ratios in solar units (Anders & Grevesse 1989) are plot- 
ted against the projected radius scaled by the virial radius, 
rigoi m all panels. 

NGC 1404 (Matsushita et al. 2007b), and a cluster of 
galaxies Abell 1060 (Sato et al. 2007a) and AWM 7 (Sato 
et al. 2008). All systems show very similar value of Si/Fe 
ratio, to be 1—1.5. Mg/Fe ratio is slightly higher in 
HCG 62, Abell 1060 and AWM 7 than in other systems, 
i.e. NGC 720, Fornax cluster, and NGC 1404. We compare 
metal abundances of NGC 507 with those in the Fornax 
cluster and HCG 62 as shown in figure 8. The solar abun- 
dance by Feldman (1992), giving [Fe/H] = 3.24 x 10~ 5 and 
employed by Matsushita et al. (2007b), was scaled by a 
factor of 0.7 to match the Anders & Grevesse (1989) value. 
The Fe abundance of NGC 507 is almost the same as that 
of HCG 62 in the central region, while the Fornax cluster 
shows a lower value. Abundance ratios of O/Fe, Mg/Fe, 



Si/Fe, and S/Fe are quite similar for these 3 systems at 
r ~ 0.1 rise- Therefore, the abundance ratios show closer 
similarity than the absolute abundance values among dif- 
ferent systems. The efficiency of the metal enrichment 
may depend on parameters such as age, starformation ef- 
ficiency, contribution from cD galaxies. However, the rela- 
tive contribution of SNe la and II and the process of metal 
mixing in the ICM seem to be quite similar for different 
clusters and groups. 

Tamura et al. (2004) reported abundance ratios for 
19 clusters studied with XMM-Newton, and the mean 
Si/Fe ratio in cool and medium temperature clusters with 
kT < 6 keV was ~ 1.4. This is consistent with the Suzaku 
results for groups and poor clusters including NGC 507. 
Their O/Fe ratio, ~ 0.6, in the cluster core also agrees 
with the Suzaku results including our NGC 507 case. 
Matsushita et al. (2003); Matsushita et al. (2007a) also 
reported abundance ratio for M87 and the Centaurus clus- 
ter, respectively, based on XMM-Newton observations. 
M87 showed Mg/O ratio to be ~ 1.3 in the central region, 
and the Centaurus cluster indicated O/Fe and Si/Fe ratios 
within 8' to be consistent with our results. The Si/Fe ratio 
for NGC 507 with XMM-Newton (Kim & Fabbiano 2004) 
and Chandra (Humphrey & Buote 2006) is ~ 1, which is 
almost the same as our result. 

Sato et al. (2007b) studied contributions of SNe la and 
II to the metal enrichment, based on Suzaku results of 
NGC 507, Abell 1060, AWM 7, and HCG 62 (Sato et al. 
2007a; Sato et al. 2008; Tokoi et al. 2008). Sato et al. 
(2007b) showed the number ratios of SNe II to la to be 
~ 3.5 for the ICM of the above systems. These results 
also suggest that the clusters and groups have passed the 
same metal enrichment process in the ICM. 

4-2. Metal Mass-to-Light Ratio 

We examined the parameter of metal mass-to-light ra- 
tios for oxygen, iron, and magnesium (OMLR, IMLR, 
and MMLR, respectively) to compare the ICM metal dis- 
tribution with the stellar mass profile. The metal mass 
profiles shown in figure 9(a) were calculated from the 3- 
dimensional gas mass profile, which is based on the surface 
brightness profile by XMM-Newton (see also figure 2 and 
subsection ), the gas density in the central region, r < 2', 
through our spectral fit as shown in table 5, and the abun- 
dance profile measured with Suzaku. The derived iron, 
oxygen, and magnesium mass within the 3-dimcnsional 
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Fig. 9. (a) Cumulative mass, M(< R), within the 3-dimcnsional radius, R, for O, Fe, and Mg in NGC 507, based on the combination 
of the abundance determination with Suzaku and the gas mass profile with XMM-Newton. (b) Ratio of the O, Fe, and Mg mass in 
units of Mq to the B band optical luminosity in units of Lq (OMLR, IMLR, and MMLR, respectively) against the 3-dimcnsional 
radius. Black symbols represent our results for NGC 507, and gray symbols of OMLR and IMLR are for the Fornax cluster by 
Matsushita et al. (2007b). 



Table 8. Comparison of IMLR, OMLR and MMLR for all systems. 





IMLR 


OMLR 


MMLR 


r (kpc/ri 8 o) 


k(T) 


Reference 


Suzaku 














NGC 720 .... 


1 x 10~ 4 


4 x 10~ 4 




25/0.04 


~ 0.56 keV 


Tawara et al. (2008) 


Fornax 


4 x 10~ 4 


2 x 10~ 3 




130/0.13 


- 1.3 keV 


Matsushita et al. (2007b) 


NGC 507 ... . 


1.7 x 10~ 3 


6.6 x 10~ 3 


1.1 x 10~ 3 


260/0.24 


- 1.5 keV 


This work 


HCG 62 


4.6 x 10~ 3 


3.8 x 10~ 2 


1.5 x 10~ 3 


230/0.21 


~ 1.5 keV 


Tokoi et al. (2008) 


A 1060 


4.0 x IO" 3 


4.3 x 10~ 2 


1.6 x 10~ 3 


380/0.25 


-3 keV 


Sato et al. (2007a) 


AWM 7 


1.1 x IO" 2 


5.8 x IO" 2 


9.4 x 10" 3 


570/0.35 


~ 3.5 keV 


Sato et al. (2008) 


XMM-Newton 
Centaurus . . . 


4 x IO" 3 


3 x IO" 2 




190/0.11 


-4 keV 


Matsushita et al. (2007a) 



radius of r < 260 kpc are 3.4 x 10 s , 1.3 x 10 9 , and 2.1 x 10 8 
Mq, respectively. 

Historically, B-band luminosity has been used for the es- 
timation of the stellar mass (Makishima et al. 2001). We 
calculated the B-band luminosity in NGC507 by transla- 
tion from the Two Micron All Sky Survey (2MASS) 2 with 
an appropriate color B — K = 4.2 for early-type galaxies 
in Lin & Mohr (2004) , along with the Galactic extinction, 
A B = 0.267, from NASA/IPAC Extragalactic Database 
(NED) in the direction of NGC 507. We used 2° x 2° data 
set centered at the NGC 507 coordinate as shown in ta- 
ble 1, and subtracted the r > 1° region as a background, 
and also removed the NGC 499 region. In addition, we 
deprojected the luminosities along with the annuli of the 
spectral fits. The resultant luminosity within this Suzaku 
observation, r < 13', is 2 x 10 11 L Q in B-band. The radial 
luminosity profile corresponds to the value of figure 9 (a) 
divided by (b). 

We calculated the integrated values of OMLR, IMLR, 
and MMLR within r < 260 kpc as shown in figure 9(b) , and 
their values turned out to be ~ 6.6 x 10~ 3 , ~ 1.7 x 10~ 3 , 
and ~ 1.1 x 10~ 3 M©/L©, respectively. The errors are 

2 The database address: http://www.ipac.caltech.edu/2mass/ 



only based on the statistical errors of metal abundance 
in the spectral fit, and the uncertainties of the gas mass 
profile and the luminosity of member galaxies are not in- 
cluded. The IMLR values are consistent with the collec- 
tive results with ASCA by Makishima et al. (2001). The 
MMLR and IMLR show similar steep increase with radius 
up to r ~ 150 kpc and seem to reach almost a plateau at 
150-300 kpc. This tendency is similar to the IMLR profile 
with XMM-Newton for M 87 and the Centaurus cluster 
by Matsushita et al. (2007a) within r < 0.1 r 180 . 

We also note that the derived OMLR and IMLR curves 
for NGC 507 in figure 9 are very similar to those for the 
Fornax cluster (Matsushita et al. 2007b) as seen in fig- 
ure 9(b). We summarize recent measurements of IMLR 
and OMLR in table 8, where one can compare the present 
results with those in Fornax cluster, Centaurus cluster 
and NGC 720. It is suggested that smaller systems with 
lower gas temperature tend to show not only lower IMLR 
in Makishima et al. (2001) but also lower OMLR and 
MMLR, even though the scatter of the data is large. 
Because O and Mg are mainly synthesized by SNe II, 
while Fe is by both SNe la and II, as shown in Sato et al. 
(2007b), OMLR and MMLR show the ratios of heavy stars 
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Fig. 10. Hardness ratio image for the two energy 
bands, 1—2 keV to 0.5-1 keV, in the central region of 
the Suzaku observation. The observed images by XISO-3 
were superposed on the sky coordinate after removing 
calibration source region, and smoothed with a Gaussian 
with a = 16 pixel ~ 17". The white box corresponds to 
the X-ray discontinuity recognized by Kraft et al. (2004) 

to galaxies. Although our observation within r < 0.3 r^o 
show no clear difference between the distribution of IMLR 
and OMLR and MMLR, we need to observe the outer re- 
gion to the virial radius in order to know the metal en- 
richment process of the ICM. 

4-- 3. Discontinuity in the Surface brightness 

Since previous X-ray observations suggested possi- 
ble dynamical features in the NGC507 system (e.g. 
citekraft04), we looked into temperature structures which 
are out of axis symmetry in the form of 2-dimensional map 
of hardness ratio. Figure 10 shows the hardness ratio im- 
age based on the intensities in 1-2 keV and 0.5-1 keV. 
Although a vignetting correction was not performed, the 
two energy bands show very similar vignetting features. 
The transmission drop due to the OBF contamination 
was not corrected for either, which is lager for the 0.5- 
1 keV band and at the central part of the CCD image. 
Therefore, the hardness ratio should show a systematic 
drop from the center to the outer regions even for a con- 
stant temperature ICM. In the hardness image, the cen- 
tral region shows clearly lower hardness, and these feature 
is reported for HCG 62 observation in Tokoi et al. (2008). 
Including this feature, deviation from the circular sym- 
metry is not significantly recognized in the hardness ratio 
image of NGC507. 

Kraft et al. (2004) found a sharp edge or discontinuity 
in the X-ray surface brightness by Chandra in the east 
and south east region at ~ 2' from the center. The dis- 
continuity corresponds to the radius where the influence of 
the central cool component becomes weak, as indicated in 



figure 10. Kraft et al. (2004) suggested three possibilities 
for the discontinuity; 1) the edge is caused by a motion 
of NGC 507 with respect to the lager scale dark matter 
potential of the group, 2) the pressure across the bound- 
ary is balanced by an unseen relic radio lobe produced 
in an earlier epoch of nuclear activity, 3) the discontinu- 
ity is created by the subsonic/transonic expansion of the 
relatively weak lobe currently observed lying interior to 
the discontinuity. Regarding the possibility 1, Kraft et al. 
(2004) noted that the temperatures inside and outside the 
discontinuity are almost the same. However, while Kraft 
et al. (2004) fitted with a single temperature model, we 
employed the two temperature model. As shown in fig- 
ure 6 and 10, relative intensity of the cool component 
drops fairly sharply with radius around the discontinuity 
region. 
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